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Abstract 

First-principles, density- functional based electronic structure calculations are carried out for 
MgC{Nii-xCox)z alloys over the concentration range < x < 1, using Korringa-Kohn-Rostoker 
coherent-potential approximation (KKR CPA) method in the atomic sphere approximation (ASA). 
The self-consistent calculations are used to study the changes as a function of x in the equation of 
state parameters, total and partial densities of states, magnetic moment and the on-site exchange in- 
teraction parameter. To study the magnetic properties as well as its volume dependence, fixed-spin 
moment calculations in conjunction with the phenomenological Landau theory are employed. The 
salient features that emerge from these calculations are (i) a concentration independent variation 
in the lattice parameter and bulk modulus at a;~0.75 with an anomaly in the variation of the pres- 
sure derivative of bulk modulus, (ii) the fixed-spin moment based corrections to the overestimated 
magnetic ground state for 0.0 < x < 0.3 alloys, making the results consistent with the experiments, 
and (iii) the possibility of multiple magnetic states at x ~0.75, which, however, requires further 
improvements in the calculations. 
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I. INTRODUCTION 

Much of the interest in the 8K cubic perovskite superconductor MgCNi^ stems from 
its spectral distribution of electronic states A sharp peak, primarily 

composed of Ni 3d bands show up in the density of states spectra, suggesting that the 
material is close to a ferromagnetic instability |2, y|- Theoretical calculations based on 
the Stoner model as well as the NMR experiments reveal the possibility of co-existence of 
spin-fluctuation with superconductivity This makes MgCNi^ a material in a class of 

its own, since according to conventional models of superconductivity magnetic scattering is 
thought to induce pair-breaking effects. A simple case study is that of the binary cubic V^A^ 
alloy in which the expected superconducting transition temperature is about 30K, which is 
lowered down to ~6K due to spin-fluctuations [fll. Ilfl|. 

The origin of the incipient magnetism in MgCNi^ could be partly attributed to the 
presence of C in its octahedral interstitial position The spatially extended C 2p orbitals 
strongly hybridize with that of the Ni 3d, de-localizing the electronic states. As a result the 
density of states at the Fermi energy, N{Ef), falls short of satisfying the Stoner criteria. 
It is suggested that approximately 0.5 holes, which may correspond to 10% of Mg replaced 
by an alkali metal, or about 7% replacement of Ni with Co, would drive the corresponding 
material magnetic 0,0,0]. The latter follows the rigid-band model, where the Fermi level 
is expected to lower in energy with respect to hole concentration, thus making Ep coincide 
with the transition-metal-derived 3d singularity. 

Previously, studies of hole dopings via vacancies and B substitutions in the C sub-lattice 

Sand Fe and Co substitutions in the Ni sub-lattice have been carried out llJ, llfiL 
0, Q, Q y , Q. I„ none of ..e above ca«e. a «„ite .a,„etic .„^t,o„ was 

found feasible. For Co substitutions the experimental results are unambiguous over the 
fact that Tc decreases as Co at% increases in the disordered MgC{Nii_xCox)3 alloys [3, 
However, they remain controversial on the rate at which the Tc is decreased 
. Resistivity characterization shows that Tc decreases gradually with increasing x. 
A small superconducting transition is observed for the samples with x=0.5, which indicates 
that the superconducting volume fraction decreases upon substitution of Co for Ni jl^ . 
However, this is inconsistent with the susceptibility measurements where Co in dilute limit 
is found to completely suppress the superconductivity, although no magnetic long-range 
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order is detected The rapid loss of superconductivity is argued to be consistent with 
magnetic quenching of the superconductivity, probably in the form of spin-fluctuations. 

Co forms a complete solid solution in MgCNi^ in the cubic perovskite structure. In 
MgC(Nii_xCox)3 alloys, the ferromagnetic onset is estimated to be around 0.75 at% of 
Co Further, experiments as well as theoretical calculations suggest that a complete 

replacement of Ni with Co drive in a definite ferromagnetic ground state |l3 . 



iJ, 3. The Co 



local moment in MqCCo^ is calculated to be 0.33/iB with the magnetization energy being 



MCf UUc 



equal to 22mRy [IJ, Theoretical calculations also reveal that the stoichiometrically 
ordered alloys MgCNi2Co and MgCNiCo2 remain non-magnetic 

As a precursor to understanding the nature of magnetic instability in MgCNi^, a the- 
oretical study of Co substitution in the Ni sub-lattice was carried out. To study the 
effects of compositional disorder on the structural, electronic and magnetic properties of 
MgC{Nii_xCox)2, alloys we have employed the Green's function technique in conjunction 



with the coherent-potential approximation |23,l2J|. The versatility of the method allows a 
relative comparison of the alloy energetics of the stoichiometrically ordered alloys and their 
disordered counterparts. 

In the present work, we have studied the changes in the equation of state parameters, 
density of states and the magnetic properties as a function of x in MgC{Nii_xCox)z alloys. 
Both self-consistent as well as the fixed-spin moment method [2^ in conjunction with the 
phenomenological Landau theory of phase transition are employed to study the magnetic 
properties of the disordered MgC{Nii_xCox)z alloys. 

The paper is organized as follows. In Sec. 2, we briefiy describe the computational 
details. Sec. 3 contains the results and discussion in terms of equation of state parameters, 
electronic structure and magnetic properties of MgC{Nii_xCox)3 alloys. 

II. COMPUTATIONAL DETAILS 

The ground state properties are calculated using the Korringa-Kohn-Rostoker (KKR) 
method formulated in the atomic-sphere approximation (ASA) [3, Q|- Chemical disorder 
in the Ni sub-lattice with Co, are assumed to be completely random and is accounted for by 
means of coherent-potential approximation (CPA) |2a|. For improving the alloy energetics, 
the ASA is corrected by the use of both the muffin-tin correction for the Madelung energy [^l 
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and the multi-pole moment correction to the Madelung potential and energy [2^, |23|. These 
corrections bring significant improvements in the calculations by taking into account the 
non-spherical part of the polarization eff'ects j^^. Note that the electrostatic interactions 
are obviously a key contribution to the eff'ective interactions between diff'erently charged 
cations. 

The partial waves in the KKR-ASA calculations are expanded up to /max =3 inside the 
atomic spheres. The multi-pole moments of the electron density have been determined up to 
^max=6, and then used for the multi-pole moment correction to the Madelung energy. The 
exchange-correlation effects are taken into consideration via the local-density approximation 
(LDA) with Perdew and Wang parametrization The core states have been recalculated 
after each iteration. The calculations are partially scalar relativistic in the sense that al- 
though the wave functions are non-relativistic, first order perturbation corrections to the 
energy eigenvalues due to the Darwin and the mass-velocity terms are included. Further, 
screening constants a and j3 were incorporated in the calculations, following the prescrip- 
tion of Ruban and Skriver j^l- These values were estimated from the order (A^) locally 
self-consistent Green's function method [^^l and were determined to be 0.83 and 1.18, re- 
spectively. The atomic sphere radii of Mg, C and Ni/Co were kept as 1.404, 0.747, and 
0.849 of the Wigner-Seitz radius, respectively. The overlap volume of the atomic spheres was 
less than 15%, which is legitimate within the accuracy of the approximation 0|. The total 
energies were calculated with 1771 k-points in the irreducible wedge of the Brillouin zone. 
The convergence in charge density was achieved so that the root-mean square of moments 
of the occupied partial density of states becomes smaller than 10~^. 

Numerical calculations of magnetic energy AE{M) for MgC{Nii_xCox)3 alloys were 
carried out at their self-consistently determined equilibrium lattice constants using the fixed- 
spin moment method of alloy theory In the fixed-spin moment method the total energy 
is obtained for a given magnetization M, i.e., by fixing the numbers of electrons with up and 
down spins. In this case, the Fermi energies in the up and down spin bands are not equal to 
each other because the equilibrium condition would not be satisfied for arbitrary M. At the 
equilibrium M two Fermi energies will coincide with each other. The total magnetic energy 
becomes minimum or maximum at this value of M. 

The diffraction experiments measure displacements as large as 0.05 A for the Ni atoms in 
MgCNi^. These, however are found to be less significant towards any qualitative analysis of 
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the materials properties. The first-principles calculations ignoring these effects, thereby as- 
suming a rigid perfect underlying cubic symmetry, find no profound ma gnet ic effects emerg- 



ma gnet ic ettects em 



ing, unlike in the case of Fe impurities in FCC Al or in BCC Zr 
Although in disordered alloys lattice relaxations around an impurity is important, for 3d 
impurities in a 3d host these effects seem to be less important 123, 1^- Usually, such 
effects stem from large atomic size mismatch of the host and the impurity atoms which is, 
however, not the case for Co and Ni ions. 



III. RESULTS AND DISCUSSION 

The MgC{Nii_xCox)z alloys are theoretically characterized in terms of the variation of 
physical properties as a function of x. These include the equation of state parameters such 
as lattice constant, bulk modulus and its pressure derivative, electronic structure expressed 
in terms of total, sub-lattice resolved partial and /— decomposed density of states. The 
magnetic properties are determined via both the self-consistent calculations and the fixed- 
spin moment calculations, including the variation in the Co local moment and the on- 
site exchange interaction constant as a function of x. For x < 0.3, the self-consistent 
calculations yield ambiguous results about the magnetic ground state of MgC{Nii_xCox)z 
alloys, therefore, we used the fixed-spin moment method in this concentration range. These 
calculations find that for low Co rich alloys, the ground state is definitely paramagnetic. In 
conjunction with the Landau theory the MgC{Nii_xCox)z alloys for 0.0 < x < 0.3 show a 
greater propensity of magnetism. However, both the self-consistent as well as the fixed-spin 
moment calculations yield consistent results for x > 0.75 alloys. 



A. Equation of state parameters 

The equilibrium lattice constants of MgC{Nii_xCox)3 alloys for < x < 1 were deter- 
mined by calculating the total energy at six lattice constants close to the expected equi- 
ibrium lattice constant and then using the third-order Birch-Murnaghan equation of state 

in 

42, |43|. Since the choice of the exchange correlation term in the Kohn-Sham Hamiltonian is 
crucial in determining the equation of state parameters we have carried out the calculations 
using three different exchange correlation functionals namely the local density approxima- 
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Figure 1: The variation in the lattice constant of MgC[Nii^xCox)z as a function of x calculated 
using the KKR-ASA-CPA method. Results for both spin polarized and unpolarized calculations 
are shown. 

tion (LDA), the generalized gradient approximation (GGA) j^l, and the local Airy gas 
(LAG) approximation j^l- One may note that the Birch- Murnaghan equation is derived 
from the theory of finite strain, by considering an elastic isotropic medium under isothermal 
compression, with the assumption that the pressure- volume relation remains linear. The 
polynomial would yield a reasonable guess to the equation of state parameters provided the 
fit is chosen in a volume range close to equilibrium. 

The equilibrium lattice constants for MgC{Nii_xCox)3 alloys decrease as x increases. 
This is consistent with the fact that the ionic radii of Co is less than that of Ni. For 
MgCNi^ the equilibrium lattice constant was calculated to be 7.139, 7.305 and 7.233 a.u.^ 
in the LDA, GGA and LAG approximations, respectively. When compared to the exper- 
iments, we find that LDA underestimates the value of the lattice constant by 1%, while 
both GGA and LAG overestimate the value. The LDA estimate matches well with the 
previous full-potential report of Singh and Mazin j^. We note that the consistency with 
the full-potential methods owes to the muffin-tin correction to the ASA, without which the 
KKR-ASA calculations find the equilibrium lattice constant for MgCNi^ as 6.986 a.u. For 
MgCCo^, the equilibrium lattice constants for both spin polarized and unpolarized calcu- 
lations turn out to be 7.078 and 7.071 a.u, respectively. A small increase in the lattice 
constant in the spin polarized calculations is consistent with the fact that magnetization 
through exchange splitting requires larger volume. 

The variation of the lattice constants as a function of x in MgC{Nii_xCox)z alloys, 
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Figure 2: The variation in the bulk modulus for MgC{Nii_xCox)3 alloys as a function of x, 
calculated using the KKR-ASA-CPA method. Results for both spin polarized and unpolarized 
calculations are shown. 

following both spin unpolarized and polarized calculations is shown in FigUl The salient 
features that emerge are (i) deviation from the Vegard's law, (ii) a weak x dependence of the 
lattice constant in the range 0.7 < x < 0.8, and (iii) deviation of lattice constants above x = 
0.7. Deviation from the Vegard's law suggests a diversity in the inter-atomic interactions. 
The Vegard's law implicitly assumes that the chemical constituents which make up the alloy 
have similar potentials and that their distribution remains truly random. Although one may 
assume Co and Ni to have similar features, the nature of interactions in them is significantly 
different because MgCNi^ is paramagnetic while MgCCo^ is ferromagnetic. This indicates 
that the C 2p hybridization with the transition metal 3d bands in MgCNi^ and MgCCo^ 
are different and that the Ni-C hybridization would be stronger than the Co-C counterpart 
in MgCCos. 

The insensitivity of the lattice constant of MgC{Nii_xCox)z alloys with respect to x in 
the range 0.7 < x < 0.8 arises due to a balance between the chemical binding forces and 
the magnetic ones. Chemical binding tends to compress the lattice while magnetization 
via exchange splitting of bands requires large volume. These two opposing forces set in to 
provide an invariance in the lattice constant for 0.7 < x < 0.8 in MgC{Nii^xCox)z alloys. 
For X > 0.7, the spin polarized calculations yield slightly higher lattice constants compared 
to the spin unpolarized calculations, which is purely due to the magnetic effects. 

The bulk modulus is proportional to the second derivative of the energy- volume curve. 
The LDA usually underestimates the lattice constant thus leading to an overestimation of 
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Figure 3: The change in the pressure derivative of bulk modulus of MgC{Nii_xCox)3 as a function 
of X calculated using the KKR-ASA-CPA method. Results for both spin polarized and unpolarized 
calculations are shown. 

bulk modulus. The emphasis here, however, is to make a qualitative judgment about the 
bulk modulus as a function of x. FiglSl shows the change in the isothermal bulk modulus of 
MgC{Nii_xCox)z alloys as a function of x for < a; < 1. We find that the bulk modulus 
increases with increasing x. 

Assuming that the alloys are isotropic and that the Debye approximation holds, the bulk 
modulus scales in proportion to the Debye temperature Qd (oc^o ) • Higher the value of 60 
higher the lattice stiffness. In MgCNi^ the frequency of a soft acoustic Ni based phonon 
mode calculated in the harmonic approximation tends to become negative 0, lill. Stabi- 
lization of this mode results in a dynamic displacement of the Ni ions perpendicular to the 
Ni-C direction which then allows each of the Ni atoms to move away from the neighboring 
C atoms towards the empty interstitial. Partial replacement of Ni with Co weakens the 
bonding with the C atoms. This reduces the advantage of relaxations or distortions thereby 
increasing B^), thus a concurrent hardening of the phonon modes. When compared to the 
spin unpolarized calculations, the spin polarized ones yield a lower bulk modulus, which is 
directly correlated with their corresponding lattice constants. 

Since the pressure derivatives of the bulk modulus involves higher order derivatives, they 
are often noisy. Qualitatively, however, in the Debye approximation the pressure derivatives 
of the bulk modulus contain the information of the averaged lattice vibrations of the material 
via the Gruneissen parameter. If one assumes that all the vibrational modes respond to 
volume in a similar fashion, then they are often useful in interpreting the structural details 
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of the material. The variation in the pressure derivative of the bulk modulus, as shown in 
FigOl suggests an anomaly at x = 0.82. The discontinuity in the case of spin polarized 
results suggests that the anomaly could be related to magnetic effects because it appears in 
the proximity of magnetic onset in the MgC{Nii^xCox)3 alloys. 





LDA (Ry) 


GGA (Ry) 


LAG (Ry) 


MgCNi2Co 


-1866.120100 


-1866.883569 


-1869.423046 


MgC{Nio.e6Coo,33)3 


-1866.118525 


-1866.882022 


-1869.421520 


AE (mRy) 


1.575 


1.547 


1.526 


MgCNiCo2 


-1815.219322 


-1815.974307 


-1818.477701 


MgC{Nio.33Coom)3 


-1815.219848 


-1815.974861 


-1818.478283 


AE (mRy) 


-0.526 


-0.554 


-0.573 



Table I: Comparison of the equilibrium total energy expressed 
in Ry units for x = | and x = ^ alloys of MgC{Nii^xCox)3 
in three different approximations to Vxc', LDA, GGA and LAG 
and their difference {AE = Eordered - Edisordered) 



The versatility of the KKR-ASA method is that it allows a direct comparison of the 
structural energies of both disordered and ordered alloys. Two ordered alloys can be identi- 
fied in the concentration range < x < 1. These are MgCNi2Co and MgCNiCo2, which 
correspond to x = 0.33 and 0.66 respectively in the case of disordered MgC{Nii_xCox)3 
alloys. Either of these alloys falls in the 0.3 < x < 0.7 range, where both spin polar- 
ized and unpolarized calculations for the disordered alloys yield degenerate total energies. 
TablelH compares the KKR-ASA structural energies for the ordered and disordered alloys, 
calculated using three different exchange correlation functionals to Kohn-Sham Hamiltonian. 
Evidently, from TableJH it follows that for x = | the ordered phases are more stable than 
their disordered counterparts. However, as x increases to |, the disordered alloys become 
more stable. The changes in the structural energies provide some clues about the chemical 
ordering of atoms due to magnetic effects. We plan to report the results in a future work. 

The change in the equation of state parameters upon chemical ordering was also calculated 
using the KKR-ASA method in conjunction with the Birch-Murnaghan equation of state 
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with LDA. The results are shown in Table|TTl The calculations shows that the disordered 
alloys have slightly higher lattice constant than their stoichiometrically ordered counterparts. 





Qeqia.U.) 


B,g{Mbar) 




MgCNi2Co 


7.113 


0.432 


4.325 


MgC{Nio,6QC 00,33)3 


7.114 


0.433 


4.326 


MgCNiCo2 


7.086 


0.451 


4.463 


MgC{Nio.33C 00,66)3 


7.088 


0.451 


4.237 



Table II: Comparison of the equation of state parameters 
for stoichiometrically ordered and disordered alloys using the 
LDA- KKR-ASA- CPA method. 



B. Electronic structure 

Electronic structure of MgC{Nii^xCox)3 alloys are characterized both by means of spin 
unpolarized and polarized density of states, calculated using the KKR- ASA-CPA method. 
For X < 0.75, spin unpolarized density of states are shown while for x > 0.75, both spin 
unpolarized and polarized density of states are shown. 

1. spin unpolarised density of states 

In FiglH we show the total and sub-lattice resolved density of states of MgCNi3 and 
MgCCo3 calculated at their equilibrium lattice constants. The isolated deep valence states 
in FiglU are essentially the C 2s states. Further higher in energy, a gap opens up followed 
by an admixture of C 2p and metal 3d states. The strong hybridization is manifested in the 
DOS as a wide energy spread of the C 2p and the transition metal 3d states. From FigHJit 
follows that the hybridization between the C-Ni orbitals are much more stronger than the 
C-Co orbitals in their respective alloys. Hybridization gives not only an important mixing 
between the states of the conduction bands but also leads to a separation of the bonding and 
anti-bonding states, creating a pseudo-gap. Existence of pseudo-gaps have been reported in 
various crystalline solids, amorphous alloysj^l and also in quasi-crystalsQ, Hi- Of the two 
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Figure 4: The spin unpolarized total and sub-lattice resolved partial density of states for MgCNi^ 
(upper panel) and MgCCo^ (lower panel) calculated at their equilibrium lattice constants. In the 
inset, a blow up around the Fermi energy for the corresponding alloys are shown. The vertical line 
through energy zero represents the alloy Fermi energy. 

mechanisms, ionic and covalent, proposed for the formation of pseudo-gap in alloys the latter 
is found to play a significant role in the case of transition-metal-based inter-metallics. The 
pseudo-gap in hexagonal close-packed transition metal rich systems are largely attributed 
to metal d resonance. Though the Mg-Ni{Co) bonding is v^^eak, a peak just below — 0.4Ry 
constitutes the bonding betv^^een 3s electrons of Mg and the transition metal 3(i electrons 
in their respective alloys. The striking difference in the DOS of MgCNi^ and MgCCo^ is 
the absence of C 2p states at the Ep for MgCCoy,. In MgCCo^ the N{Ep) is entirely due 
to the Co 3d states vi^hile for MgCNi^ the C 2p contribution to N{Ep) is about 7%. 

The present KKR- ASA-CPA calculations find N{Ep) for MgCNi^ to be 14.56 state/Ry 
atom. The value is, hov^^ever, at variance v^^ith the previous reports. For example, Szajek 
reports the value as 14.32 state/Ry atom j^, Mazin and Singh report as 13.57 states/Ry 
atom y|. Shim et al as 14.52 states/Ry atom Ijl, Rosner et al as 13.06 states/Ry atom 

Q R n 

Dugdale and Jarlborg p| report N{Ep) to be 17.27 and 9.49 states/Ry atom p| for 
two different band-structure methods with exchange-correlation effects considered in the 
LDA and using the experimental lattice constant. These results show that N{Ep) is indeed 
sensitive to the type of the electronic structure method employed and also to the numerical 
values of the parameters like that of the Wigner-Seitz radii and others. Note that these 
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differences are significant as tfiey control tfie proximity to magnetism in tfie Stoner model, 
as emphasized by Singh and Mazin For MgCCo^ the value of N{Ep) is determined to 
be 14.43 state/Ry atom, in agreement with the earlier reports. 

FigEl shows the spin unpolarized, total and sub-lattice resolved partial DOS for 
MgC{Nii_xCox)3 alloys for x= 0.2, 0.4, 0.6 and 0.8. The overall characteristic features 
of the DOS structure are preserved as x increases in MgC{Nii^xCox)3, however with the 
sharp structures being smeared out due to disorder. As x increases, the decreasing electron 
count is followed by the inward movement of Ep on the energy scale. The bottom of the 
states, however, remain more or less fixed with respect to the alloy Ep as a function of x. 
At about X = 0.7 the Ep is pinned in the pseudo-gap. The resulting N{Ep) is minimum for 
these alloys over the whole concentration profile. As x increases further from 0.7, the N{Ep) 
gradually increases due to the Co 3d states. A high N{Ep), predominantly metal 3d in char- 
acter, thus increases the possibility of exchange splitting of bands in MgC{Nii_xCox)3. 




Energy (Ry) 

Figure 5: The spin unpolarized total (solid line), sub-lattice resolved C (shaded) and {Nii-xCox)3 
density of states, with x as shown in the panels, for MgC{Nii-xCox)3 alloys. The vertical line 
through energy zero corresponds to the alloy Fermi energy. 

One may note that the C 2p states in and around Ep decrease with increasing x, and 
virtually cease to exist for x > 0.7 alloys. Since the number of C electrons remains invariant, 
one may then naturally expect a significant charge redistribution on the energy scale. This 
is illustrated in FigEl Though the overall shape of the DOS remains the same, a significant 
charge redistribution over the energy range —0.6 Ry to —0.2 Ry is observed. 
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Figure 6: The spin unpolarized sub-lattice resolved C 2p DOS for MgC{Nii-xCox)3 with x as 
shown in the panel. The vertical line through energy zero corresponds to the alloy Fermi energy. 
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Figure 7: The variation in the total density of states at Fermi energy as a function of x in 
MgC{Nii^xCox)3 for < x < 1 alloys in units of states/Ry cell. 

in FigUl we show the evolution of the states at Ep as a function of x in MgC{Nii^xCox)z 
alloys. Initially with increase in x, the N{Ep) slightly increases, however, not as expected 
from the rigid band model of alloy theory. Disorder smears out the Ni 3d- derived van Hove 
peak of MgCNi^ considerably. 

Thus the essential difference between the DOS of MgCNi^ and MgCCo-s is the absence of 
C contribution to N{Ef). One can find that for both MgCNi^ and MgCCo^ the N{Ef) are 
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comparable and are equal to 14.557 states/Ry atom and 14.425 states/Ry atom, respectively. 
However, their properties are found to be significantly diff'erent. MgCNi^ is non-magnetic 
with incipient magnetism akin to spin-fluctuations being anticipated, while MgCCo^ is a 
non superconductor with a definite ferromagnetic state. 

The change in the total density of states for the chemically ordered MgCNi2Co and 
MgCNiCo2 alloys with respect to their disordered counterparts, MgC^NioseCoo,^^)^ and 
MgC{Ni(),33Coo.m)3 respectively, are shown in FiglHl The overall characteristic features 
of the DOS remain essentially the same, with exception to the peak broadening in the 
disordered alloys. The change in the Ni 3d and Co 3d states for the corresponding ordered 
and disordered alloys are shown in FiglHl and the change in the C partial DOS is shown in 
FigHni 
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Figure 8: A comparison of the spin unpolarized total DOS of stoichiometrically ordered and disor- 
dered MgC{Nii-xCox)3 alloys for x= | and |, respectively. The vertical line represents the Fermi 
energy. 



2. Spin polarized density of states 

With N{Ef) in the paramagnetic DOS showing high values for x > 0.75 and with no C 2p 
states at or near Ep, a possibility of electronic distortion leading to a ferromagnetic ground 
state can well be anticipated. Taking the cue, spin polarized calculations were carried out 
for MgC{Nii_xCox)3 alloys for 0.7 < x < 1.0. 
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Figure 9: A comparison of spin unpolarized, sub-lattice resolved transition metal DOS of stoichio- 
metrically ordered and disordered M gC {Nii^xC Ox)z alloys for x= ^ (upper panel) and | (lower 
panel). The vertical line represents the Fermi energy. 

Fig. [TT] shows the spin polarized, total and sub-lattice resolved partial DOS of MgCCo^. 
The deep states are primarily C 2s states. Higher in energy and close to Ep, the states are 
predominately Co 3d in character. As expected, the Co 3d bands are exchange split due to 
magnetic effects. As a measure of the exchange splitting on the energy scale, we find the 
difference between the centres of the 3d up bands and that of the down bands C^ as 22 
mRy, which corresponds to a local magnetic moment of 0.33 fiB- These values are consistent 
with earlier full-potential calculations. Small exchange splitting of the Mg and C 2p bands 
are seen due to hybridization effects, however, the splittings are negligible when compared 
to that of the Co 3d bands. 

The change in the DOS as a function of x in MgC{Nii^xCox)z alloys are shown in FigIT2l 
For both majority and minority bands, the deep states essentially consist of C 2s states. 
Higher in energy towards Ep-, a gap opens up followed by an admixture of C 2p - metal 
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Figure 10: A comparison of spin unpolarized sub-lattice resolved partial C 2p DOS of stoichio- 
metrically ordered and disordered M gC {Nii^xC Ox)z alloys for x= | (upper panel) and | (lower 
panel). The vertical line represents the Fermi energy. 
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Figure 11: The spin polarized total and sub-lattice resolved partial spin up and down DOS of 
MgCCo^. The vertical line through the energy zero represents the Fermi energy and + and - DOS 
represent the majority and minority spin DOS respectively. 



3d states. The spin up bands are fully occupied and the Ep for these states are pinned in 
the pseudo-gap. For the spin down states, the Ep moves slightly inward with increasing x. 
Magnetism is thus governed by the change in the position of Ep in the minority DOS with 
respect to that of the majority bands. 

The sub-lattice resolved, concentration weighted, /—decomposed 3d DOS of 
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Figure 12: The spin polarized total density of states of M gC (Nii^xC Ox)3 alloys for 0.85 < x < 
0.95. The vertical line through the energy zero represents the Fermi energy and + and - DOS 
represent the majority and minority spin DOS respectively. 

MgC{Nii_xCox)z alloys are as shown in FigUHl where the exchange splitting as well as 
the distortions can be seen clearly. It also becomes evident that not all d sub-states are 
equally exchange split. Relatively, dy^ and (is^a.i bands show the maximum exchange split- 
ting, while the states contained in the x — y plane, the distortions in the majority and 
minority bands are more noticeable. Distortions result due to significant charge transfer. 
Bonding is anisotropic due to covalent character along the x — y plane and rather metallic 
character along the z direction. 



C. Magnetic properties 

1. Variation of Co local moment 

In general, when a paramagnetic DOS shows a high N{Ep), with Ep in the anti-bonding 
region, the system is close to a magnetic instability. Under these conditions, distortion in- 
evitably decreases the total energy. The first choice, however, is a reshuffling of the electronic 
states instead of the atoms getting displaced. The electrons prefer to rearrange themselves 
resulting in a spontaneous magnetization thereby lowering the electronic symmetry and 
annihilating the anti-bonding states. 
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Figure 13: The spin polarized, sub-lattice resolved and orbital decomposed Co 2>d partial density 
of states for 0.85 < x < 0.95 alloys. The vertical line through the energy zero represents the Fermi 
energy. 

According to Dronskowski et al 0|, anti-ferromagnetism is likely to step in when 
the paramagnetic DOS has Ep in a non-bonding region. For those systems, where Ep is in 
the anti-bonding region, exchange splitting of bands towards ferromagnetism is more likely. 
In MgCNi^ the Ep is in the anti-bonding region, and hence it is closer to ferromagnetic 
instability rather than an anti-ferromagnetic one. Further, anti-ferromagnetism had been 
ruled out in MgCNi^ owing to the absence of nesting features in the Fermi surface topology. 

The present calculations find small moments at Co site for x < 0.3, as shown in Fig. [TH 
The results are inconsistent with experiments. Some of the factors that can influence the 
magnetic moment at the Co site are the choice of the sphere radii, inadequate basis set used 
to expand the wave functions, lack of lattice relaxation and local environment effects. In this 
direction, calculations were carried out for different sphere radii, but the local moment at 
the Co site remained inevitable. To include local environment effects, super-cell calculations 
based on the LMTO-ASA as well as locally self consistent Green's function method were 
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employed, however the Co remained magnetic. Since common to all the above calculations 
was LDA, the appearance of local magnetic moment at Co site was thought to be due to 
the limitations of LDA itself. 




0.4 0.6 

Co at% 



Figure 14: The variation in the Co local magnetic moment in M gC {Nii^xC Ox)3 alloys calculated 
using the KKR-ASA- CPA method in LDA for < a; < 1. 



It is known that first-principles LDA based calculations for materials such as FeAl j^, 
NisGa H, Sr3Ru20j 



SrRhO^ and NaQ,^Co02 yield a finite local moment 

in disagreement with the experiments which find them to be paramagnetic. The error that 
evolved due to LDA in the present self-consistent calculations for MgC{Nii^xCox)'3 alloys, 
was however corrected by means of the fixed-spin moment method, the results of which we 
show below. 



2. On-site magnetic exchange interactions 

The total exchange coupling parameter Jq is calculated as a function of x, x > 0.7 

n 

following the prescription of Liechtenstein et al [60]. Essentially the method is based on 
mapping the change in energy due to the deviation of a single spin at a reference site "o", 
from the collinear ferromagnetic ground state onto an effective Heisenberg model. The 
change in the energy corresponding to this small spin density perturbation is approximated 
by the change in the sum of one electron energies by appealing to Andersen's local force 
theorem. Using Lloyd's formula one can express the sum of one electron energies in terms 
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of scattering path operator or the auxiliary Green's function. The latter can be evaluated 
from the knowledge of potential functions and the structure constant matrix. Mapping the 
change in energy onto a Heisenberg model results in an expression for the exchange coupling 
constant: 

i 

which is the sum of exchange interactions between the reference spin and all its neighbors. 
In a mean field theory this coupling constant is proportional to the Curie temperature Tc of 
the system. For a multi-component random alloy the mean field estimate of the Tc can be 
assumed to be the concentration-weighted average of the coupling constants calculated for 
the component atoms. Thus, for the case oi MgC{Nii^r^Cox)z alloys, the average coupling 
constant can be calculated as 

Jo = (1 - x) Jm + {x) Jco 
where Ji {i = Ni or Co) is given as 

= E P dE Im {Al(^) [Tl^iz) - Tl^{z)] + Al{z)Tl^{z)A^^{z)Tl^{z)} 

with 

- PUiz) - Pl^iz) 
PUz) ^[z- CQ [Ai, + 7i. {z - CI,)] 

TLi^) = {9Laiz)) {l + (PLaiz) - PlM ){gL.{z))} 

where a is the spin index (t or J,), and P^o- ^^e potential function of the component i for 
the orbital L and spin a. The potential function has been expressed above in terms of C, 
A and 7 of the LMTO Hamiltonian, z is the complex energy and Ep is the Fermi energy of 
the alloy, (giaiz)) is the configurationally averaged auxiliary Green's function within CPA 
and PLa{z) is the coherent-potential of the medium. The CPA calculations are performed 
by invoking site-site approximation. 
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Figure 15: The variation of the on-site exchange coupling constant Jq in niRy, as a function of 
X > 0.75 in MgC{Nii^j,Cox)3 alloys. 

FigflBl shows the variation in the Jm and Jco as a function of x in MgC{Nii_xCox)z 
alloys. Both Jco and ^co (ref. Fig lT^ follov^ a similar trend. Both shov^^ a slov^ variation in 
and around the ferromagnetic onset followed by a sharp increase, further getting saturated 
when x— >1. Except for an overestimation of the exchange coupling constant, the calculations 
are expected to produce the correct trend. The disappearance of the average Co local 
moment is accompanied by disappearance of the average exchange interaction J^o, a result 
which can be related to the strong magneto- volume effect. 



3. Fixed-spin moment calculations 

The density functional theory is an exact ground state theory. However, common ap- 
proximations like LDA to the density functional theory is thought to miss certain important 
physics due to the assumption of a uniform electron gas especially when spin-fluctuation 
effects are important. The effects of spin-fluctuations are, however, described on a phe- 
nomenological level using the Ginzburg-Landau theory. Although a robust quantitative 
theory based on first-principles is yet to be implemented, it is possible to make an estimate 
of spin-fiuctuation effects based on LDA fixed-spin moment calculations. 

According to Mazin and others 0|, the overestimation of the tendency of metals to- 
wards ferromagnetism within the LDA can be used as an indicator of critical fiuctuations 



21 



in a material. However, for this to be an effective indicator, competing states, like anti- 
ferromagnetism need to be ruled out in each material. As mentioned above, since Ep resides 
in the anti-bonding region as well as the absence of nesting features in the Fermi surface 
topology of MgCNi^ shows that the material under consideration may be far from any 
anti-ferromagnetic instability. 

For MgC{Nii_xCox)z, the overestimation of the moments for x < 0.3 could be taken as 
an indication of spin-fluctuations, following the works of Singh and Mazin j^. However, no 
attempt is made to determine the magnitude of the fluctuations in the present work. The 
study is limited to understanding the propensity of magnetism with Co substitutions, for 
which one can use the phenomenological Ginzburg- Landau coefficients, the input to which 
follows from the fixed-spin moment method. 

Numerical calculations of magnetic energy Ai?(M) for MgC{Nii_xCox)z alloys in the 
range < x < 1 are carried out by the fixed-spin moment method. One can then write 
down a Ginzburg-Landau expansion for the magnetic energy, Ai?(M) = I]">i ^a2jM^' for 
n = 3. The variation of Ai?(M) with respect to M is shown in Fig^j For x < 0.3, the 
curves are relatively fiat near M = and fit well to the form as given above with n = 3. 
The sensitivity of A£'(M) with M becomes more distinguishable as x increases from 0.3 
and beyond, i.e., the fiatness of the curve disappears. For x > 0.7 the magnetic energy is 
negative, indicating a lower energy ferromagnetic state to be stable. 

Fig [121 shows the variation of Ginzburg-Landau coefficients 02 and 04 as a function of 
X in MgC{Nii_xCox)3 alloys. The coefficient 02 is of prime importance and is used as 
a precursor to indicate the magnetic phase of the system. A positive 02 represents the 
paramagnetic state while a negative 02 indicates the ferromagnetic state. From Fig IT7I we 
find that initially, as x increases, 02 tends to zero, but remains positive. This indicates that 
a ferromagnetic state is unstable for low Co rich alloys. However for x > 0.75, 02 is negative 
indicating a ferromagnetic state for the system. These calculated results are now consistent 
with the experiments, which show a definite paramagnetic phase for low Co rich alloys and 
a magnetic state for high Co rich alloys. 

As the unit cell volume is increased, Ep is lowered in energy. This will then sharpen the 
structures in the DOS. If N{Ef) then satisfies the Stoner criterion, i.e., I*N{Ep) > 1, where 
/ is the Stoner exchange integral, then the system would become magnetic. Having seen 
that the self-consistent calculations yield ambiguous magnetic properties, one can extend 
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Figure 16: The variation in the AE{M) as a function of magnetization M, calculated at the 
respective equilibrium volume for MgC{Nii-xCox)3 alloys as a function of x, as indicated in the 
figure. 
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Figure 17: The variation of Ginzburg-Landau coefficients as a function of Co concentration (x) in 
MgC{Nii_xCo^)3 alloys. 

the fixed-spin moment calculations to the magnetization dependence of MgC{Nii_xCox)z 
with respect to the unit cell volume. 

Fig^l shows the variation of the magnetic energy as a function of volume, expressed in 
terms of its equilibrium volume Vq, for x = 0.01,0.1 and 0.2, respectively. For x = 0.01, an 
expansion by 8% shows no magnetic phase for the alloy. However, for both x = 0.1 as well 
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Figure 18: The variation of AE{M) as a function of M for MgC{Nii^xCox)z alloys at various 
volume ratios as indicated in the figure for (a) x = 0.01 (upper panel), (b) x = 0.1 (middle panel) 
and (c) X = 0.2 (lower panel). 

as X = 0.2, a magnetic phase looks possible just above its equilibrium volume. The curves 
dip into the negative region of the magnetic energy as can be seen from the middle and the 
bottom panel of Fig. ^1 respectively, where the ferromagnetic state is lower in energy. 

For a; = 0.1 and 0.2, a ferromagnetic phase becomes stable just above its equilibrium 
volume suggesting a magneto- volume instability in MgC{Nii^,j.Cox)z alloys. Fitting the 
above data to the Ginzburg-Landau expression for n = 3 for the three concentrations x = 
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0.01, 0.10 and 0.20, one finds that for expanded volumes a magnetic phase becomes definitely 
possible as Co at% increases in MgC{Nii_xCox)3. 
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Figure 19: The variation of the Ginzburg-Landau coefficients, a2(in units of a4(in units of 



^), and oe (in units of ^) as a function of volume for MgC{Nii-xCox)3 alloys for (a) x = 0.01 
(upper panel), (b) x = 0.1 (middle panel) and (c) x = 0.2 (lower panel). 

FigHni shows the variation of the Ginzburg- Landau coefficients as a function of volume 
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ratio for x — 0.01, 0.10 and 0.20, respectively. For ^ slightly greater than unity, for both 
X — 0.1 and x — 0.2, the coefficient 02 reverses its sign, indicating a phase transition with 
respect to volume. This confirms that these alloys are on the verge of a magneto-volume 
instability. The higher order coefficients play a significant role in determining the alloy 
magnetic properties. 

It has also been noted that the accuracy of the total energy, evaluated by means of ASA 
is limited to a few mRy. From calculations, one finds that the change in the energy is far too 
small for any reliable quantitative analysis. However, within the approximations mentioned 
and with well defined k— mesh in the Brillouin zone, one can certainly infer that the trend 
produced via the fixed-spin moment method for MgC{Nii_xCox)z alloys are qualitatively 
correct. 

IV. CONCLUSIONS 

For low Co rich MgC{Nii-xCox)^ alloys, we find the possibility of enhanced spin- 
fiuctuations in the material, which is close to a ferromagnetic instability. However, for 
Co rich alloys, a definite ferromagnetic ground state exists. At the magnetic cross-over, 
i.e., for X = 0.75, a concentration independent variation in the structural properties are 
determined. For example, the lattice constant and the bulk modulus appear to be constant 
under the given approximations of the KKR- ASA-CPA theory. The pressure derivative of 
the bulk modulus indicates, within the Debye approximation, a significant change in the 
averaged vibrational modes. Our calculations show that the electronic structure of the dis- 
ordered MgC{Nii_xCox)z alloys deviates significantly from that of the rigid band model. 
The striking feature is the recession of C 2p states towards lower energies as a function 
of increasing x in MgC{Nii_xCox)z alloys. The self-consistent calculations overestimate 
the magnetic moments at the Co site for low Co rich alloys, due to the limitations in the 
local-density approximation. Corrections to the overestimation in the magnetic moments 
are accomplished by means of fixed-spin moment method, in conjunction with the Ginzburg- 
Landau free energy functional. It then follows from the fixed-spin moment calculations that 
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for X < 0.7 the material is definitely paramagnetic. 
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